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Abstract 

In  GPS  one  of  the  primary  errors  contributing  to  positioning  inaccuracy  is  the  performance  of 
the  on-board  atomic  clock.  To  determine  and  predict  the  performance  of  this  atomic  clock  has 
been  a  problem  due  to  the  ambiguity  of  the  orbital  position  error  and  clock  uncertainity  in  the 
Radio  Frequency  (RF)  tracking  of  the  navigation  signals.  The  Laser  Retrorefiector  Experiment 
(LRE)  on-board  NAVSTAR  35  and  36  provides  a  means  of  separating  these  ambiguious  errors 
by  enabling  highly  precise  and  accurate  satellite  positions  to  be  determined  independent!)'  of  the 
RF  signals.  The  results  of  examining  onboard  clock  behavior  after  removing  the  orbital  position 
signatures  mil  be  discussed.  GPU  RF  tracking  data  from  various  DoD  and  other  sites  are  used  to 
reconstruct  the  onboard  clock  data  and  examine  the  clock  behavior.  From  these  data,  the  effects  of 
clock  performance  on  GPS  positioning  performance  can  examined. 


INTRODUCTION 

The  purpose  of  this  project  is  to  identify  and  investigate  means  of  enhancing  the  Global 
Positioning  System  (GPS)  system  integrity  and  performance.  This  project  involves  installing  laser 
retrorefiector  array.s  onDboard  Global  Positioning  System  (GPS)  satellites,  tracking  the  satellites 
involved  in  cooperation  with  the  NASA  Satellite  Laser  Ranging  (SLR)  network  and  collecting 
these  data  for  analysis  and  comparison  with  GPS  pseudorange  data.  The  Laser  Retrorefiector 
Experiment  (LRE),  previously  known  as  the  Advanced  Clock  Ranging  E.xperiment  (ACRE)91, 
was  .submitted  by  the  U.  S.  Naval  Research  Laboratory  (NRL)  to  the  TriDService  Space  Test 
Program  for  spacecraft  integration  funding  as  a  triDservice  space  experiment.  The  objective 
of  such  an  experiment  is  to  provide  an  independent  high  precision  measurement  to  compare 
or  calibrate  the  GPS  pseudoDranging  signal.  This  project  is  a  cooperative  effort  involving 
the  NASA  Goddard  Spaceflight  Center  SLR  group,  the  NRL  and  the  University  of  Maryland. 
Installation  of  the  LRE  on  the  GPS  .satellite  was  performed  in  conjunction  with  the  GPS  Joint 
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The  GPS  system  is  a  predicted,  realDtime,  passive  ranging  navigation  system,  made  up  of  space, 
control  and  user  segments.  The  space  and  control  elements  comprise  the  system  proper,  and 
the  user  segment  operates  passively  utilizing  the  products  of  the  system  transmitted  by  the  space 
segment.  The  user’s  information  is  computed  from  the  control  segment’s  tracking  network’s 
data  and  other  data  provided  by  external  sources,  such  as  the  U.S.  Naval  Observatory  (USNO) 
for  Universal  Coordinated  Time  (UTC)  corrections.  The  tracking  network  data  are  similar  in 
content  to  that  used  by  the  user  segment  and  is  relayed  to  the  Master  Cxrntrol  Station  (MCS) 
for  computation  and  prediction  of  the  .system  states  which  are  uploaded  into  the  satellites  for 
the  users.  Embedded  in  the  space  and  control  segments  are  atomic  clocks  to  maintain  all 
elements  of  the  system  in  synchronization.  These  atomic  clocks  enable  the  precise  time  of 
propagation  measurements  (known  as  Pseudoranges)  the  users  measure  to  determine  range 
between  themselves  and  the  satellites,  and  the  capability  of  determining  the  precise  positions  of 
the  satellites  needed  as  the  users’  position  reference.  Small,  passive  LRE  on  two  GPS  satellites, 
capable  of  supporting  highly  precise  laser  ranging  to  that  satellite,  tracked  by  a  worldwide 
network  of  SLR  stations  are  to  produce  highly  precise  and  accurate  orbital  ephemerides.  These 
data  are  being  compared  with  GPS  orbits  generated  by  the  MCS  and  the  Defense  Mapping 
Agency  postDprocessed  precise  ephemerides  to  separate  the  satellite  position  and  onDboard 
atomic  clock  errors.  This  error  separation  should  provide  a  foundation  for  better  understanding 
the  satellite  clock  onDorbit  performance,  error  propagation  within  the  MCS  data  computation 
process,  and  an  independent  calibration  of  GPS  accuracy. 

SATELLITE  EQUIPMENT 

The  LRE  is  a  panel  of  a  laser  retro  reflector  cubes,  24  x  19.4  cm  (9.45  x  7.64  inches)  as 
shown  in  Figure  1.  This  array  consists  of  32,  2.7  cm  (1.06  inch)  reflectors  of  the  design  used 
onDboard  Glonass  satellites.  These  arrays  were  built  and  tested  by  the  Russian  Instutite  for 
Space  Device  Engineering  in  a  cooperative  arrangement  with  the  University  of  Maryland.  The 
placement  on  the  selected  satellites,  NAVSTAR  35  and  36,  is  shown  in  Figures  2  and  3. 

LASER  TRACKING  NETWORK 

The  laser  returns  from  the  LRE  is  estimated  to  be  a  factor  of  36  lower  than  that  of  Glonass, 
whose  array  size  is  about  120  x  120  cm  (47.2  x  47.2  inches),  and  a  factor  of  3  to  4  lower 
than  Etalon  (the  Russian  laser  retroreflector  satellite  at  Glonass/GPS  altitudes).  Good  Glonass 
returns  to  the  NASA  mobile  laser  sites  (MOBLAS)  are  roughly  equal  to  that  from  LAGEOS. 
LAGEOS  is  routinely  tracked  by  the  NASA  and  cooperating  laser  sites.  For  Etalon  tracking,  a 
receiver  threshold  of  4  photo  electrons  is  used  by  MOBLAS  for  day/night  operation.  With  the 
LRE  and  the  same  receiver  threshold,  the  ranging  returns  are  estimated  to  be  10  to  20return, 
ranging  returns  could  be  increased  to  about  the  same  level  as  Etalon  if  the  receiver  thresholds 
on  the  MOBLAS  were  reduced  from  four  photoelectrons  to  one  photoelectron  (lunar  mode) 
during  nightDtime  tracking.  Daylight  tracking  from  MOBLAS  is  more  difficult  due  to  the  high 
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background  noise  rate  and  the  single  stop  time  interval  units  used  rather  than  the  multistop 
event  timers  used  at  the  lunar  ranging  sites.  Modifications  to  enable  daytime  tracking  from 
MOBLAS  has  been  prototyped  and  proven  at  the  GSFC  tracking  .site  and  the  MOBLAS  .sites 
are  being  upgraded. 

The  results  presented  here  are  for  NAVSTAR  35  only.  NAVSTAR  36  was  launched  significantly 
later  and  has  only  been  sporatically  tracked.  There  are  twelve  sites  which  with  varied  frequency 
have  successfully  tracked  NAVSTAR  35.  The  U.S.  systems  at  Monument  Pk.,  CA,  Greenbeit, 
MD,  Quincy,  CA,  McDonald  Obs.,  TX,  Haleakala,  HI,  Yarragadee,  Australia  and  the  inter¬ 
national  sites  at  Herstmonceux,  U.K.,  Graz,  Austria,  Wettzell,  Germany,  Potsdam,  Germany, 
Maidanak,  Uzbekistan  and  Evpatoria,  Ukraine.  The  distribution  of  the  tracked  “segments”  by 
each  of  these  stations  indicate  that  some  of  the  sites  have  only  tracked  over  certain  periods  of 
time  in  a  non-uniform  way.  This  is  due  to  the  fact  that  tracking  has  been  limited  to  daylight. 
Consequently,  there  are  only  short  periods  of  a  day  or  so  when  several  sites  were  simultaneously 
successful  in  tracking  the  satellite.  In  particular,  on  November  18,  1993  ten  passes  of  data 
were  acquired.  This  is  the  rea.son  why  this  day  was  chosen  to  do  preliminary  comparisons  with 
the  GPS-derived  orbits  for  NAVSTAR  35. 

GPS  TRACKING 

For  intercomparison  with  the  GPS  derived  data,  these  data  are  being  collected  at  NRL  along 
with  the  laser  tracking  data.  Tracking  data  from  the  GPS  Control  Segment  stations,  USNO, 
the  broadcast  position  data  and  DMA  precise  ephemerides  are  being  collected.  These  data  are 
continuous  over  the  inDorbit  operation  of  the  satellites.  To  utilize  the  GPS  derived  tracking 
data  for  intercomparison  with  the  laser  derived  data,  the  local  clocks  at  the  GPS  Monitor 
Station  sites  must  be  accounted  for  since  they  are  the  basis  for  the  GPS  tracking  measurements. 
In  GPS  itself  these  clocks  are  accounted  for  by  the  use  of  GPS  Time  which  is  a  common 
synchronization  time  computed  at  the  MCS.  However,  the  GPS  ranging  measurements  are 
directly  related  to  the  local  clocks  whose  performance  must  be  removed  if  the  satellite  clock  is 
to  be  isolated  from  the  satellite  orbital  position  and  evaluated.  The  laser  data  is  independent 
of  this  influence  on  ranging  measurements  since  the  local  clock  i.s  used  for  timetagging. 

To  determine  the  performance  of  the  station  clocks,  common  view  time  comparisons  with 
USNO  were  made  to  the  Colorado  Springs,  Hawaii  and  Ascension  stations.  These  comparisons 
provide  local  station  clock  compared  to  the  Master  clock  at  USNO.  These  data  show  that  large 
jumps  and  discontinuities  are  present  as  shown  in  Figures  4,  5,  and  6.  These  jumps  are  due  to 
changes  in  the  local  docks  or  switching  necessary  for  the  operation  of  the  system.  Navigation 
users  would  not  be  aware  of  these  changes  since  they  use  GPS  Time  which  is  a  computed  time 
accounting  for  these  changes.  For  this  experiment,  removal  of  the  local  clock  and  the  satellite 
position  error  by  laser  derived  positions  from  the  GPS  tracking  data  will  leave  the  satellite 
clock  as  the  principal  error  component. 
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ORBITAL  ANALYSIS  AND  RESULTS 


The  lERS  Standards!^*  with  minor  excursions  (e,g.  JGM-2  gravity  field  vs.  GEM-T3)  have 
been  adopted  to  ensure  as  much  compatibility  with  other  analyses  results  as  possible.  The  orbits 
are  integrated  in  the  mean  system  of  J2000  and  only  the  terrestrial  effects  due  to  relativity 
are  used.  Modeling  of  the  perturbing  forces  on  the  satellite  is  tailored  after  the  LAGEOS 
SLR  analysis  standards.  The  exception  is  the  limited  gravity  field  terms  (18,18)  required  here 
due  to  the  higher  orbit  of  the  target  satellite.  The  time-varying  part  of  the  geopotential  is 
accommodated  by  modeling  the  solid  Earth  and  oceans  tidal  accelerations  and  the  secular 
change  in  the  terrestrial  oblateness.  Because  the  NAVSTAR  satellites  are  not  passive  as 
LAGEOS,  attitude  variations  must  be  accounted  for  and  the  implications  these  have  on  the 
solar  and  thermal  forces  acting  on  the  satellite  a  various  times.  The  model  used  to  describe 
these  forces  is  the  abridged  version  of  Rockwell  International’s  ’'ROCK42’'  model  by  Fliegel- 
Gallini-Swift,  the  T20f51.  An  additional  acceleration  along  the  satellite  body-fixed  Y-axis,  the 
so-called  Y-bias,  is  also  adjusted.  Due  to  the  length  of  the  arc  used,  once  per  revolution 
accelerations  (with  constraints)  are  also  included  and  adjusted  over  the  same  intervals  that  the 
constant  accelerations  apply.  The  duration  of  these  intervals  is  variable  and  they  have  been 
kept  constant  as  long  as  the  data  allow  in  order  to  increase  the  robustness  of  the  .solution.  The 
strategy  followed  has  been  to  keep  the  same  number  of  adjusted  accelerations  while  lengthening 
the  arc  and  to  introduce  a  new  set  of  accelerations  once  the  data  indicate  a  change  in  the  orbit. 
These  parameters  along  with  the  state  vector  at  epoch  are  the  only  force  model  parameters 
that  are  adjusted. 

Measurement  modeling  accounts  for  tropospheric  refraction,  tidal  variations  of  the  site  including 
ocean  loading  (in  all  three  directions),  tectonic  motions,  and  occasionally  measurement  biases. 
The  tropospheric  refraction  model  for  SLR  is  the  Marini-Murray  model.  Ocean  loading 
effects  at  the  SLR  sites  was  computed  using  the  Scherneck  model  for  the  eleven  main  tidal 
constituents  of  Schwiderskii’s  ocean  tidal  model.  Tectonic  motions  for  the  sites  are  either  from 
the  LAGEOS-based  solution  SL8.3f*l  or  the  NDVEL-lNNRl‘^1,  Only  simple  measurement 
biases  were  adjusted  on  a  few  occasions  for  certain  sites.  Most  of  these  biases  are  the  result 
of  “fine-timing”  of  the  ranging  gates  at  the  site  in  order  to  achieve  the  maximum  number  of 
returns  possible.  Once  the  sites  are  equipped  with  the  better  detection  packages  there  should 
be  no  need  to  change  these  thresholds  and  therefore  the  chance  of  introducing  biases  to  the 
data  will  be  minimized. 

The  collected  SLR  data  are  analyzed  and  reduced  based  on  the  force  and  measurement  models 
described  in  the  previous  section.  A  long  arc  of  about  104  days  was  continuously  extended  as 
new  data  become  available.  This  arc  was  used  to  check  on  the  fidelity  of  the  force  model. 
The  data  fit  the  arc  with  an  rms  of  3  cm.  The  geographical  distribution  of  the  data  set  did  not 
include  southern  hemisphere  tracking  and  that  can  introduce  significant  biases  in  the  orbits. 

Table  1  shows  the  rms  residual  for  each  of  the  tracking  sites.  It  is  hard  to  assess  the  quality  of 
the  orbits  without  a  uniform  data  distribution.  November  18,  1993  being  the  best  tracking  day 
within  our  data  set,  it  was  used  as  a  test  day  to  verify  orbit  quality  and  gain  some  insight  in 
the  level  of  agreement  with  the  “‘radiometric  data”  — determined  orbits  that  the  International 
GPS  Service  (IGS)  for  Geodynamics  is  routinely  distributing!*!.  Two  fourteen  day  arcs  were 
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fit  to  the  data;  one  for  November  5-18  inciusive  and  one  beginning  on  November  18.  These 
arcs  have  only  12  hours  worth  of  data  in  common;  11:00  UT  to  23:00  UT,  on  November  18. 
The  data  fit  either  arc  with  an  rms  residual  of  about  1.9  cm.  In  both  cases,  the  state  vector 
and  one  set  of  accelerations  were  estimated.  The  two  orbits  are  based  on  just  over  200  normal 
points  each.  For  arcs  of  such  length  this  can  hardly  be  called  a  sufficient  amount  of  data.  The 
trajectories  from  the  two  adjustments  were  then  compared  in  terms  of  radial,  cross-track,  and 
along-track  differences  over  their  common  segment.  The  statistics  from  this  comparison  (mean 
and  rms  about  the  mean),  are  shown  in  Table  2. 


Table  1 

Residual  statistics  for  the  104-day  SLR-determined  arc 


Site 

No.  of  Obs. 

RMS  (cm] 

Monument  Peak,  CA 

311 

2.3 

Haleakala,  HI 

215 

3.1 

McDonald  Ob.s.,  TX 

81 

2.7 

Quincy,  CA 

4 

0.1 

Greenbelt,  MD 

8 

1.0 

Graz,  Austria 

175 

2.8 

Herstmonceu.x,  U.K. 

101 

3.4 

Potsdam,  FRG 

47 

2.1 

Wettzell,  FRG 

121 

3.1 

Totals 

1063 

2.9 

Table  2 

Trajectory  Differences  for  the  two  SLR-determined  14— day  arcs. 

Component 

Position  [cm] 

Velocity  [cm/s] 

Cross 

Cross  I  Along 

Mean 

5.1 

21.8 

-19.0 

0.0028 

0.0002  0.0012 

RMS 

3.2 

37.0 

10.9 

0.0017 

0.0015  1  0.0059 

Despite  the  fact  that  the  SLR  data  distribution  is  not  as  optimal  as  would  be  preferred  for  a 
precise  orbit  determination,  it  is  still  worthwhile  comparing  to  the  GPS-derived  orbits  distributed 
by  IGS  for  geodetic  work.  The  IGS  orbit  was  rotated  into  the  inertial  frame  and  used  as 
‘'observations”  with  the  GEODYN  data  analysis  software  package  to  restitute  a  dynamic  orbit 
fitting  that  data.  The  converged  trajectory  was  then  compared  to  the  SLR-derived  orbit  in  the 
radial,  cross-track,  and  along-track  directions  (Figure  7).  Statistics  of  these  differences  of  the 
IGS  orbit  from  both  SLR  14-day  arcs  are  shown  in  Tables  3  and  4.  The  common  segment  of 
course  is  only  one  day  (November  18)  in  both  cases. 


Table  3 

Trajectory  Differences  SLR— 1 

vs.  IGS  GPS  orbit 

Component 

Position  [cm] 

Velocity  [cm/s] 

Direction 

Radial 

Cross 

Along 

Radial 

Mean 

8.9 

RMS 

7.7 

EEB 
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Table  4 

Trajectory  Differences  SLR-2 

vs.  IGS  GPS  orbit 

Component 

Position  [cm] 

Velocity  [cm/s] 

Direction 

Cross 

Along 

Radial 

Cross 

Along 

Mean 

3.6 

41.5 

58.7 

RMS 

9.8 

72.9 

CONCLUSIONS 

The  collection  of  the  GPS  tracking  data  is  proceeding  well  and  the  SLR  data  is  proceeding 
slowly.  The  complication  of  removing  the  local  atomic  clock  offset  and  drift  from  the  GPS 
data  is  being  accomplished  using  the  common  view  technique  of  simultaneous  observations  of 
the  satellites  at  two  sites.  These  comparisons  should  be  of  sufficient  accuracy  to  remove  these 
effects  from  the  individual  satellite  tracking  data.  With  SLR  derived  positions  having  sufficient 
confidence  the  resulting  satellite  atomic  clock  performance  should  be  isolated  for  evaluation. 

With  limited  SLR  data,  it  is  hard  to  come  to  firm  conclusions.  The  two  orbit  comparisons 
show  at  least  the  level  of  compatibility  of  the  SLR  and  IGS  orbits  at  about  10  cm  in  the  radial 
direction,  whether  it  be  in  the  mean  or  the  rms  sense.  This  is  a  very  limited  test,  where  neither 
technology  has  put  forward  its  best  accomplishments  and  capabilities.  A  much  more  uniform 
and  extended  SLR  data  set  will  be  required  before  we  can  reliably  determine  an  orbit  at  the 
few  centimeter  level  of  accuracy.  On  the  other  hand,  reduction  of  GPS  data  directly  within 
GEODYN  will  remove  any  inconsistencies  in  the  standards  and  the  reference  frame  used  by  the 
IGS  analysis  centers  and  the  SLR  group.  Upcoming  modifications  to  the  SLR  ground  receivers 
will  allow  for  a  further  increase  in  the  tracking  capabilities  of  several  additional  sites  and  add 
the  needed  southern  hemisphere  tracking.  An  initial  effort  to  compare  the  SLR  derived  orbits 
with  those  distributed  by  IGS  indicates  that  the  two  agree  at  the  decimeter  level  radially  and  at 
the  0.5-1. 0  meter  level  in  the  cross-track  and  along-track  directions.  The  amount  of  collected 
data  by  site  and  geographical  region  is  far  from  optimal  for  a  reliable  orbit  determination,  so 
these  results  should  be  interpreted  with  caution. 
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QUESTIONS  AND  ANSWERS 


MARC  WEISS  (NIST):  On  one  of  the  plots  of  the  residuals,  I  wasn’t  exactly  sure  what  the 
data  meant.  There  were  normal  plots  for  the  laser  ranging,  and  I  thought  they  were  open 
squares.  Were  those  DMA  or  isiso-ephemeris  ranging? 

RONALD  BEARD  (NRL):  The  normal  points  from  the  satellite  data  you  mean? 

MARC  A.  WEISS  (NIST):  Yeah. 

RONALD  BEARD  (NRL):  I  think,  as  John  mentioned  yesterday,  they  are  doing  a  number 
of  pulses,  like  10  pulses  per  second,  to  get  the  returns.  They  have  taken  like  five  minutes  of 
these  returns,  and  they  averaged  those  into  one,  what  they  call  a  “normal  point.” 

MARC  A.  WEISS  (NIST):  And  you  were  comparing  those  on  the  same  plot? 

RONALD  BEARD  (NRL):  The  normal  points  are  made  to  the  raw  range  measurements,  if 
that  is  the  one  I  think  you  mean. 

MARC  A.  WEISS  (NIST):  It’s  the  first  one.  And  then  there  was  an  RMS  of  some  two 
millimeters.  The  open  squares  are  what? 

RONALD  BEARD  (NRL);  The  open  squares  are  the  raw  range  measurements  that  they 
are  making.  They  are  getting  like  10  a  minute,  or  10  a  second. 

MARC  A.  WEISS  (NIST):  So  the  RMS  is  really  the  self-  consistency  of  the  range  mea¬ 
surements  with  the  laser. 

RONALD  BEARD  (NRL):  That’s  correct. 

MARC  A.  WEISS  (NIST):  Okay.  I  understand  that  you’re  trying  to  do  orbit  reconstruction 
based  on  laser  measurements  only.  And  it  seems  that  you  can  get  a  simple  measure  of  the 
consistency  by  just  looking  at  range  measurements  for  your  laser  and  range  estimates  from, 
say,  DMA  orbits  or  broadcast  orbits.  Has  that  been  done? 

RONALD  BEARD  (NRL):  Yes  and  no. 

MARC  A.  WEISS  (NIST):  That  seems  a  lot  simpler,  I  would  be  very  interested  to  know 
how  they  compare  simply  for  range  measurements. 

RONALD  BEARD  (NRL):  It’s  a  lot  more  difficult  than  it  appears  on  the  surface.  That’s 
one  of  the  reasons  we  want  to  try  to  do  some  simultaneous  tracking,  so  we  can  do  just  that. 
Even  the  locations  of  the  stations  and  the  lasers,  it’s  difficult  to  get  enough  correlation  between 
the  two  to  just  simply  do  a  comparison  of  those  two.  But  we  have  been  trying. 

JOHN  LUCK  (ORRORAL  OBSERVATORY):  First  remark:  I  think  the  comparison 
between  the  SLR-derived  orbits  and  the  IGS  orbits  for  35  and  36  are  consistent  at  about  15 
to  20  cm  level.  The  graph  that  you  were  just  looking  at  is  the  self-consistent  residuals  for  the 
laser-derived  orbit. 

My  question  was:  Seeing  that  this  is  a  very  powerful  tool  for  geodetic  investigation,  such  as 
height  determinations,  sea-level  monitoring  and  things  like  that,  are  there  any  plans  to  include 
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retro-reflector  arrays  on  future  GPS  spacecraft?  And  if  so,  could  you  please  make  them  bigger? 

RONALD  BEARD  (NRL):  Well,  no  and  yes.  There  are  no  plans  to  include  them  downstream 
that  Pm  aware  of.  There  are  no  specific  plans.  There  are  recommendations  for  doing  that, 
and  various  options  have  been  discussed.  If  we  do,  we  sure  have  the  world  as  our  incubator. 
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